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Covariant Constitutive Relations, Landau Damping and Non-stationary 
Inhomogeneous Plasmas 
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UK 

Models of covariant linear electromagnetic constitutive relations are formulated that have wide applicability to 
the computation of susceptibility tensors for dispersive and inhomogeneous media. A pcrturbative framework 
is used to derive a linear constitutive relation for a globally neutral plasma enabling one to describe in this 
context a generalized Landau damping mechanism for non-stationary inhomogeneous plasma states. 
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Constitutive relations are widely used when describing 
the behaviour of electromagnetic fields in continuous me- 
dia. Although their domain of applicability is often deter- 
mined experimentally, causality and locality play impor- 
tant roles in their theoretical foundations. Limitations 
arise due to the inherent non-linearity contained in the 
classical equations describing the coupling between the 
motion of individual particles or continuous charge dis- 
tributions and a self-consistent electromagnetic field. For 
small disturbances perturbative linearisation techniques 
are available and approximation schemes exist for calcu- 
lating effective susceptibility tensors that arise from these 
constitutive relations. The effective constitutive relations 
that result from such schemes often rely for their valid- 
ity on assumptions such as material homogeneity and 
non-relativistic perturbations about stationary configu- 
rations. This letter addresses some of the issues that arise 
when some of these assumptions are relaxed and the de- 
gree to which concepts such as Landau damping can be 
generalized for relativistic inhomogeneous plasmasi^^. 

The macroscopic Maxwell equations can be written as 
the exterior system 



dF = 0, 



e d * G = - * J 



(1) 



where F is the Maxwell 2-form, G = eo F + II the excita- 
tion 2-form and J the source current vector on a space- 
time M, in terms of the Hodge map * and metric dual 
J associated with the spacetime metric g. In this form 
certain 4-current sources — d * II arc included in G with 
the remaining 4-currents denoted by J. A covariant con- 
stitutive model provides functional relations between G 
(or II) and F and between J and F. Relative to any 
unit future-pointing timelike observer vector field U on 
M, the forms F and G define the frame dependent elec- 
tromagnetic 1-forms e u = ijj F, b u — ijj * F, d u = ijjG 



and h u = i v * G so that F = U A e u - *(U A b u ) and 
likewise for G. 

If one restricts to causal linear responses a natural co- 
variant constitutive relation is given by the non-local ex- 
pression 



J A 



U[F] ab (x) = 



y<£j-{x) 



Xa b c d (x,y)F ed (y)dy cde f (2) 



where Xabcd(x,y) is a two-point susceptibility kernel. 
The events x and y are given in arbitrary coordinates 
with summation over Latin indices from to 3 and 

dydef _ ^yC ^ ^yd ^ ^ye ^ ^y f _ rp ne q^jJ structure 

has been imposed by requiring that Xabcd(x,y) = if y 
does not lie in the past light cone, J~(x), of x. This con- 
stitutive relation can be used to model media which are 
spatially inhomogeneous and temporally non-stationary 
and is meaningful in spacetimes containing gravitation. 

To facilitate the discussion and the use of two-point 
tensors introduce two copies of M, denoted Mx and My, 
with generic points x G Mx and y € My coordinated 
by (x°,...,x 3 ) and (y°,...,y 3 ) respectively. The val- 
ues Xabcd(x,y) denote the coordinate components of the 
4— form field x over the product manifold Mx x My in 



the induced coordinates (a; , . . . , x , y , . . . , y ): 

X = jXabcd dx a A dx b A dy c A dy d . 



(3) 



In terms of x an d the projection py : Mx x My — > My, 
py(x,y) — y equation ([2]) can be written 



U[F] = [ 

J M Y 



XApyF. 



(4) 



The tensor x has 36 independent components since dx ab 
and dy cd are antisymmetric so Xabcd = ~Xbacd = ~Xabdc- 
A special case of |2]) arises in Minkowski spacetime. 
Being parallelizable it admits a family of translation maps 
A z : M —> M, x i-)- A z (x) = x + z for all events x in M. 
This induces the translation maps B z : Mx x My — » 
Mx x My, B z (x,y) = (x + z,y + z). Imposing this 
transnational symmetry on %, i.e. B*x = Xi equation ([2]) 
can be written 



Ul[F]ab(x) 



yeJ-(x) 



X abc d{x - y)F ef {y)dy 



cdef 



(5) 



where X abcd (x~y) = Xabcd(x,y) = Xabcd{x + z , y + z) for 
any z. Since IItL4*-F] = AJIIifF] this relation describes 
spatially homogeneous and temporally stationary media 
and remains non-local in both space and time. Such 
a medium exhibits both spatial and temporal dispersion 
as follows: Minkowski spacetime admits preferred global 
Lorentzian coordinate systems, (x 3 ), with as- 

sociated cobases (dx° , dx 1 , dx 2 , dx 3 ) in which the compo- 
nents of the metric are diag(— 1, 1, 1, 1). If one defines for 
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any scalar <f> the Fourier transform 

£(fc) = / e~ ihx <j>{x) dx 0123 

with respect to such a coordinate system then 

V4F}ab( k ) = X abcd (k)F ef (k)e cde f in terms of the stan- 
dard constant alternating symbol e cde f . Such a relation 
can give rise to dispersion in media. 

If the bulk 4-velocity of a non-accelerating medium is 
V, with constant components in the above coordinate 
system, i.e. W — 0, a particular model for X in ([5|) is 
given by 

Il Ia [F} ab (x) = f V(x-y)(i v FAV) ab (y)dy 0123 

-*x [ M(x-y)(i v * Y FAV) ab (y)dy 0123 

(6) 

where V and M. are polarization and magnetization sus- 
ceptibility scalars respectively. The Fourier transform 
of ^ yields the simple constitutive relations for a spa- 
tially and temporally dispersive homogeneous isotropic 

medium: d v a {k) = (e + V(k))e%{k) and h v a {k) = 

(/i^ 1 + M(k))b^(k). If V is not inertial (W ^ 0) then 
([6]) is not a special case of ([5]) and its Fourier transform, 
although local in k, is not of the form above. 

For media that lack spatial dispersion the history of 
the medium may give rise to temporal dispersion alone. 
This can be expressed geometrically in terms of tensor 
transport along the integral curves C x : K — > M of the 
4-velocity field V of the medium. If these curves are each 
parameterized by proper time r let ^(x) be a map that 
transports tensors at C x (t) to tensors at C x (f) along 
each integral curve of V. Natural choices of transport 
maps include Lie, parallel (with respect to some space- 
time connection V) and Fermi-Walker transport. Dif- 
ferent choices of ^(x) correspond to different electro- 
magnetic responses of the medium to the disposition of 
the integral curves of V in the spacetime history of the 
medium. If Y(z) denotes a tensor field mapping 2-forms 
at z to 2-forms at z, a constitutive relation for a spatially 
inhomogeneous medium may be written 

Th[F\(C a (r)) = [ T Y*(T,T, X )($Ux)(F(C x (t))))dT 

J — CO 

(7) 

where F*(t,t,x) = (x)(Y(C w (t - f))) is a tensor 
at C x (t). This is another special case of ©. Since 

n n [F*](c x (f)) = K(*)(n*[F}{c x (T))) 

where F*(C w (t)) = ^ t {x){F{C x {t))) the medium is 
said to be stationary with respect to the transport map 
and hence V and ([7]) is valid in any spacetime. This gen- 
eralizes the notion of a temporally stationary medium in 



a spacetime with timelike Killing vectors. The tempo- 
ral dispersive properties of the medium are best defined 
with respect to a modified Fourier transform that remains 
valid in a general spacetime and employs the transport 
map along the family of curves describing the history of 
the medium. For any tensor field a and curve C x define, 
at the event C x (0), the tensor: 

/oo 
e- iuT $° T (x)(a(C x (T)))dT. 
-oo 

Then then constitutive relation II]i[F](w, x) — 
Y(cd,x)fF(u),x)) describes an anisotropic, spatially 
inhomogeneous but temporally dispersive medium. 

If V is geodesic (i.e. Vy V = 0) and &f(x) describes 
parallel transport then a particular model for Y in is 
given by 

Tlj La [F](C x (T)) = 

J T V{T~T,x){^l(x){i v FAV){C x {T))yr 

-* f M(T-T,x)($$(x)(i v *FAV)(C m (T)j\dT 

where V and A4 are spatially inhomogeneous polariza- 
tion and magnetization susceptibility scalars respectively. 
This describes non-magneto-electric, spatially inhomo- 
geneous but temporally dispersive media with consti- 
tutive relations d v (uj,x) — (eo + V(uj, x))e v (uj, x) and 

W(u,x) = {jiQ 1 +M{uj,x))W{u,x). 

The historic covariant constitutive relations proposed 
by Minkowski belong to a class of local linear relations 
on spacetime of the form 

n m [F](x) = Z(x)(F(x)) (9) 

where the tensor Z (x) maps 2-forms at events x to 2- 
forms at x. For a medium with bulk 4-velocity field V this 
tensor may be chosen so that d v (x) = (eo + V(x))e v (x) 
and h \x) = (McT 1 + M(x))b v (x) valid for inhomoge- 
neous and non-stationary isotropic media in an arbi- 
trary spacetime. However in a Minkowski background 
their Fourier transforms for non-constant V and M. yield 
non-local constitutive relations among the Fourier com- 
ponents of the electromagnetic fields so do not describe 
normal dispersive continua. 

The different constitutive relations ©, ©, ©, © 
and ([5]) above are applicable to the phenomenological 
description of linear media that exhibit temporal and/or 
spatial dispersion but rely on some knowledge of the elec- 
tromagnetic response of systems in either inertial or co- 
moving reference frames. However such relations do not 
encompass the effective susceptibility that arises when 
one applies a perturbative analysis to processes involv- 
ing non-stationary, spatially inhomogeneous plasmas^—. 
Such processes are not uncommon in astrophysical ap- 
plications or in regimes where instabilities arise^r— from 
inhomogeneities such as those in laser-plasma systems. 
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In these situations the microscopic Maxwell system for a 
neutral plasma, dF = 0, e$d * F = — * J, can be cou- 
pled with equations for the charged sources J in order 
to effect a linearisation. Since d * J — it is convenient 
to write *J = —d * II and eo F = G — II, to cast the 
system into the form used above for neutral polarizable 
media. Writing II as a functional of F in the form ([2]). 
enables one to derive an effective constitutive relation for 
a globally neutral plasma. 



I. THE COVARIANT SUSCEPTIBILITY KERNEL FOR A 
NON-STATIONARY, SPATIALLY INHOMOGENEOUS 
PLASMA 

Consider a macroscopically neutral plasma composed 
of particles labeled by the species index |pf|, mass m la] 
and charge q lc "\ described dynamically by the coupled 
relativistic (collisionless) Maxwell- Vlasov equations for F 
and the one-particle probability distributions / Lcvl (x, v) in 
an arbitrary (background) gravitational field: 



(10) 



dF = and e d*F=-*J (11) 

with the total current 1-form J(x) = g a bJ b (x) dx a given 
by 

j\x) = -J2l lal J V " lde * 9]1/2 f l "\x,v)dv™ (12) 



and Liouville vector field W lal 
d 



(13) 



— + (- T\ f { X )v°vf + l^F ef (x) 9 -vf) — 



with summation over Greek indices from 1 to 3. The 
function v°(x, v 1 , v 2 , v 3 ) is a solution of v a v b g a b{x) = — 1. 

A linear constitutive relation arises from a perturba- 
tion of this system about a background (zeroth-order) 
spatially inhomogeneous and temporally non-stationary 
solution, Fq(x), /q" 1 (x, v). The standard perturbation ex- 
pansion about such a solution yields a linear system for 
fi and Fi in terms of /o and Fq that can be solved in 
principle by the method of characteristics"^. This yields 
a first order system of integro-differential equations for 
Ff. 



dFi = 0, d-k Fi = -d-klL^Fi] 



(14) 



where LTi is a linear functional depending on /o and the 
solutions to the zeroth-order Lorentz force equations 



q 



m 



(15) 



after elimination of f±. This gives rise to a class of solu- 
tions xi to (H]) where II, x are replaced by Hi,Xi respec- 
tively. 

If one employs standard inertial coordinates in a 
Minkowski spacetime and a zeroth-order electromagnetic 
field Fq — 0, in zeroth-order all particles move along 
straight line time-like geodesies independent of |a] and 



x (2dx 0fl + e d,Tjk e^ U(T u v Uddx jk ) A dy lh 



where f(x,y) — (~g(x — y, x — y)) 1 ^ 2 and u(x, y) — 
[x — y)/f(x, y) which is manifestly not a function of x — y 
alone. 



II. LANGMUIR MODES FOR A NEUTRAL 
NON-STATIONARY, SPATIALLY INHOMOGENEOUS 
PLASMA IN MINKOWSKI SPACETIME 

The general solution to the second equation in (|T4| is 
given by eo*Fi = — -kILi[Fi] +d(3 where j3 is an arbitrary 
1-form. We define the generalized Langmuir sector to 
contain particular solutions satisfying 



eoFi = -n x [Fi]. 



(16) 



This then reduces to the standard perturbative solution 
describing longitudinal plasma oscillations about a sta- 
tionary /o in a spatially homogeneous plasma. More 
generally consider the situation where Fq = and pla- 
nar inhomogeneities arise from the following zeroth-order 
non-stationary spatially inhomogeneous solution to the 
Maxwell- Vlasov system: (|104I13[) 

(17) 



If 



m Lio£|(l + u 2)l/S 



where ,4L ion l(£) normalizes CE7|), then /L ion l initially at 
t = represents a distribution of ions where, at each 
spatial point £, the velocities belong to the 1-dimensional 
Maxwell- Jiittner distribution, but where the temperature 
TL lon l(£) and the number density of ions nL lon l(£), de- 
pend on position. It follows from (fl7|) that also ini- 
tially represents a position dependent Maxwell- Jiittner 
distribution where n^(£) = n^^) and T^(£) = 
T^(^)m^/m^ on \ After the initial moment, the ions 
and electrons drift according to (IT7)) and velocities do 
not remain in the Maxwell- Jiittner distributions. Alter- 
natively (|17p might describe a plasma composed of par- 
ticles and anti-particles. 
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In the theory of a homogeneous plasma one has the 
solution F\ = dt A d£ E(lo, n) e~ lut+1K ^ provided to and 
k satisfy a transcendental dispersion relation. This re- 
lation contains an integral that is potentially singular. 
The Landau prescription circumvents this singularity by 
complexifying lo and defining an analytic continuation 
for the integral in the complex lo plane. In an inhomo- 
geneous plasma there is no such time harmonic solution 
or associated algebraic dispersion relation between lo and 
k. We therefore propose solving (IT51) with a longitudinal 
field Fi represented as 



F x (t, = dt A d£ 



duj 



u)— — oo 



dK E(u, K)e- %u)t+i ^ 



In this case the Landau dispersion relation is replaced 
by an integral equation for E(lo, k). This equation con- 
tains a double integral that requires analytic continua- 
tion in the complex u) plane for its definition. If one 
restricts to modes with real n there are now two sin- 
gular branch points in the lo plane at lo = db|/c|. To 
define an analytic continuation from the upper-half lo 
plane the associated cuts are disposed along the half- 
lines {u = —\k\ — iX, A > 0} and {lo — \k\ — iX, A > 0} 
(see figure). 
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The integral equation for E(lo, k) is now given by 
Qo f°° 

E(lo,k) = / P(L0 7 K,k)dk 

if Im(w) > or |Re(w)| > \k\ , 
P(lo, k, k) dk — iQo / R(u),K,k)dk 



E{lo,k) 



Qo 

2tt 



Z7T 



if Im(cj) < and Re(w) < \k\ 

P(LO,K,k)dk / R(ui,K,k)dk 

00 ^ J —00 

if Im(cj) = and |Re(w)| < \k\ 



where Qo = 9 Ao^ 1 ' 31 , the principal value integral 



is given by 

P{lo, k, k) — 

~ kl E(u+k',k)(k-k)' 



-\ K -H\ (wk -u>k- k'k) 2 
and the residue by 



h K — k, 



(( K -k) 2 -K> 2 ) 



1/2 



dK' 



R(lo, k, k) 



\k-k\ BE . 



K\K\ 0lo\ K 'Tr K ' ( K 2_ W 2)1/ 2 J 

( wk A dh / _„ s K s K ^ k Lo \ 



(k 2 - w 2 ) 3 / 2 



Here s K = k/\k\ and h(n,u) — f°l e lKS h(s,u)ds. 

The square root (k 2 — w 2 ) 1 / 2 is defined so that for lo 6 K, 
\lo\ < \k\ and with the branch cuts given in the figure then 
(k 2 -lo 2 ) 1 / 2 >0. 

For real k these integral equations can be analyzed 
numerically 12 in the different domains in the lo plane and 
wave instability is associated with solutions for which 
Im(u;) > 0. Although the nature of Landau damping 
(Im(w) < 0) in the presence of inhomogeneities is clearly 
more complicated that analogous damping in homoge- 
neous plasmas, the results above indicate how the mech- 
anism depends on the nature of the initial state and an- 
alytic continuation in the complex lo plane. 
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